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Abetract 

A finite element velocity potential program 
has been developed for NASA Lewis nt the Georgia 
Institute of Technology to study acoustic wove 
propagation in complex geometries. For irrotatlon- 
al flows, relatively low sound frequencies, and 
plane wave input, the finite clement solutions show 
significant effects of inlet curvature and flow 
gradients on the attenuation of a given acoustic 
liner in a realistic variable area turbofan inlet. 
In addition, as shown in the paper, the velocity 
potential approach can not be used to estimate the 
effects of rotational flow on acoustic propagation 
since the potential acoustic disturbances propagate 
at the speed of the media in sheared flow. Ap- 
proaches are discussed chat are being considered 
for extending the finite element solution to in- 
clude the far field as well as the internal portion 
of the duct. A new matrix partitioning approach is 
presented that can be incorporated in previously 
developed programs to allow the finite clement cal- 
culation to be marched into the for field, The 
partitioning approach provides a large reduction in 
computer storage and running times. 


Introduction 

When using numerical techniques in calculating 
acoustic waves propagating in ducts (such ns finite 
elements in Refs. 1 to 3 and finite differences in 
Refs. 4 to 6) the velocity potential formulation of 
the acoustic wave equation* offers many advantages 
over the conventional linearized gas equation. z 
For two dimensional flows, the velocity potential 
approach reducea the computer storage and running 
time by an order of magnitude compared to the more 
general linearized gas equation approach. Since 
the flow into on inlet is usually modeled by poten- 
tial flow (excluding the boundary layer), the acous- 
tic velocity potential is Ideally suited for acous- 
tic Inlet calculations. 

On the other hand, when rotational flow exists 
in the inlet (wall and centerbody boundary layers), 
the accuracy of a potential flow calculation is 
questionable. The first purpose of this paper was 
to investigate the limitations of shear on the po- 
tential Clow formulation and to suggest a means of 
extending the potential flow finite element analy- 
sis to incorporate sheared flow. A second purpose 
of this paper was to determine the sensitivity of 
an acoustic liner to variations in area and flow 
experienced in a typical engine inlet.? This was 
accomplished using the finite element velocity po- 
tential computer program developed by the Georgia 
Institute of Technology for NASA Lewis Research 
Center written under NASA Lewis grant (NAS-3036). 

Finally, finite elements can simultaneously 
include both the far field and internal ducting of 
a turbofan engine, For example, in the absence of 
flow, a formulation of a finite element approach 


for sound propagation in the near field of an in- 
finite baffle with an embedded circular piston was 
presented in Ref, 8. Although the finite element 
approach is relatively easy to formulate, core 
storage requirements and long run times presently 
prevent the application of the finite element ap- 
proach to simultaneously Including botli the inlet 
and for field in a single calculation. In Ref. I, 
for example, by structuring the column matrix of 
modal unknowns to successively be the real and then 
the Imaginary part of the acoustic potential at a 
node, the global matrix is well bonded. Even so, 
the formulation of Ref. 1 was limited to 270 ele- 
ments, Reference 3 Improved the matrix structure 
by using a complex solution technique which reduced 
the bandwidth of the matrix oven further; however, 
150 total elements were the maximum quoted in the 
paper. 

The third purpose of this paper, therefore, is 
to discuss current work for extending finite ele- 
ment velocity potentials from the Internal portion 
of the duct into the far field. In particular, a 
new numerical partitioning method is presented 
which is suited for extending the solution from the 
Internal portion of the duct into the far field. 


Nomenclature 

C], constant of integration, N/m z (Eq. (22)) 

C 2 >C 3 constants of integration, m/sec (Eqs, 23 
and 24) 

c 0 velocity of sound, m/sec 

D 0 diameter at rotor position, m 

f frequency, Hz 

It height of duct, m 

t /n: 

K number of grid points In axial direction 

L length of duct, m 

M Mach number of axial mean flow, U/c 0 

m spinning mode number 

p acoustic pressure, p(x,y,t), N/m z 

p 0 entrance acoustic pressure, p 0 (o,y,t), 

N/m z 

p° spatial acouatlc pressure, p°(x,y), N/m z 

r radial coordinate, m 

t time, sec 

U axial mean flow velocity, U(x,y), m/sec 

u acoustic axial velocity, u(x,y,t), m/sec 

u° spatial acoustic axial velocity, u°(x,y), 

m/sec 

V transverse mean flow velocity, V(x,y), 

m/sec 

v acoustic transverse velocity, v(x,y,t), 

m/sec 

v° spatial acoustic transverse velocity, 

V°(x,y) , m/sec 
x axial coordinate, m 

y transverse coordinate, m 

Z impedance, g/cm z sec 

Z e ff effective impedance, g/cm z sec 

z axial position 

^ specific acoustic impedance 


r| dimensionless frequency D 0 f/c 0 or llf/c 0 

0 angular position, radians 

p 0 density, g/m 3 

‘1 steady potential function, $(x,y), m^/soc 

ft* potential function, ^(x.y.t) , m^/sce 

9 acoustic potential function, 9(Kiy»t)> 

m 3 / sec 

9° apetlal acoustic potential function, 

(f°(x,y), m^/sec 

aj angular frequency, rad/sec 

Subscripts; 

( ) x , 1 derivatives with respect to x (similar 
( )xx] for y and t) 


Sheared Mean Flow 


Rotational velocity fields exist for sheared 
viscous flows In the Inlet along the walls of a 
turbofan cowl nnd along Its conterbody. The com- 
plete acoustic equations for sheared viscous flows 
were developed by Hunger nnd Gladwell. 3 By common 
practice, the viscous terms in the linearized acous- 
tic equotlona arc neglected and the lnvlscld acous- 
tic equations arc solved using the mean flow veloc- 
ities that result from a consideration of viscosity 
only on the mean flow field, Similarly, because of 
the significant computational advantage of the 
acoustic potentiol formulation, It would be desir- 
able if the existing potential wave program could 
be extended to include boundary layer flows. First, 
the fully Irrotatlonal acoustic flow equations will 
be presented. Next, this paper will develop the 
acoUBtic flow equations for a rotational mean flow 
field with an irrotatlonal acoustic field, The con- 
sequence of such an assumption is discussed and 
conclusions are made about the use of the acoustic 
potential function ta model sheared flows. 


Irrotatlon a l Mean and irrotatlonal Acouatlc 
Flow Fields 

For two dimensional flows, die velocity poten- 
tial representation of the lnvlscld momentum nnd 
mass continuity equations can bo written aa (Ref, 10, 
pg. 76). 


<lct + 2i’ x '!$t + 20yi£y + « x + 

+ - c2<.* x + c|4.* y (1) 

where 9* 1 b the velocity potential. The symbols 
are defined in the list of symbols. To keep the 
analysis simple, the rectangular coordinate system 
Is used rather than a more general coordinate sys- 
tem. The velocity potential l 1 * is assumed to be 
composed of a steady potential 4 1 and an acoustic 
potential tp: 

4>*(x,y,t) = ^(x.y) +9(x,y,t) <2) 

The mean flow velocities are of the form 

U = 0 X = U(x,y) (3) 

V = 4>y = V (x,y) (4) 


tlon, may be regarded os negligible compared to the 
remaining terms, Eq. (1) reduces to 

(1 - M2)<P XX + <P yy - 2 ^ <P xC 
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+ 2 »i ’xy + ^g’y < 5 > 

c o 

where the mean flow velocities are related by the 
Irrotatianallty condition 


dM dV 
° dy dx 


( 6 ) 


Parallel shear flow with M(y) a nd V ° 0 is 
used to approximate boundary Layer flows. In this 
case, Eq. (5) reduces to 


(1 - n2)<p xx + q) yy - Z Ji <p xt 

' 2M § v y * "2 ’tc - 0 < 7 > 

c o 

However, Eq. (7) violates the irrototionality as- 
sumption which was used to establish Eq . (1) since 
Eq, (6) is no longer satisfied. Could, however, 

Eq. (7) be used to approximate a shear flow field? 
Equation (7) is similar to the Pridmore-Brown 
(Ref. 11) shear flow equation, except the Pridmore- 
Brr-n equation has two dependent variables. When 
the wove equation is expressed in one dependent 
variable, It becomes third order (Ref. 12, Eq. 1.22, 
pg. 9). Therefore Eq. (7) may alter the true phys- 
ics of the acoustic propagation. It would he dif- 
ficult to evaluate the effect of this approximation 
from numerical solutions. Consequently, a deriva- 
tion will now be made assuming an acoustic poten- 
tial In a parallel sheared mean flow so that Eq. (7) 
can be evaluated. 


Rotational Mean Flow and Irrotatlonal 
Acoustic Propagation 

Because of the significant computational ad- 
vantages of the acoustic potential formulation, an 
attempt will be made to extend the existing poten- 
tial wave program to include the case of boundary 
layer flows thut exist along the engine walls and 
centerbodics. To do this In a rigorous manner, one 
starts with the linearized gas dynamic equations 
for a rotational flow field. 

For the simple case of a rectangular coordi- 
nate system and isentropic (p = c^p) parallel 
sheared mean flow (1) " ll(y); V « 0), the govering 
gas dynamic equations 33 become 

Continuity 




and a constant speed of sound c Q is assumed. 

With these assumptions, and by dropping those terms 
which according to the assumption of small perturba- 
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A solution to Eqs. (8) to (10) will, of course, 
require greater computer storage and operating time 
than Eq. (7) since die number of dependent vari- 
ables has increased from one (y) to three (p,u,v). 
Now, die mean flow is considered to be rotational 
and the acoustic field is assumed to bo irrotafion- 
nl ouch that 


^2 
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This assumption reduces Eqs, (8) to (10) to 


(U) 

(’ 2 ) 


ft + U + c o p o^xx + %) - 0 


4’tX + WPx* + <Py j~ - - ^ 5^ 


l An 


<Pty + ^xy = -^f 


(13) 

(14) 

(15) 


Equations (13) to (15) can not be combined to 
obtain a second order differential wave equation in 
the dependent variable. However, differentiating 
Eq , (14) with respect to y and Eq. (15) with re- 
apect to x and aubtrncting one equation from the 
otlior yields , 

CP** + 'PyP ff + <Py $ - 0 < 16 > 


which is the governing equation for describing an 
irrotar.ional acoustic disturbance ir o parallel 
mean shear flow. Since Eq. (16) contains no deriva- 
tives with respect to time, the irrotationality as- 
sumption changes the nature of the propagation 
equation. 


The significance of Eq. (16) can be clearly 
understood by considering the simple case of a 
linear shear profile (c)2uA)y2 = 0), In this olin- 
pler but practical case, Eq. (161 reduces to 


'Pxx + 'Pyy “ 0 

or if a harmonic time dependence is assumed 
rp = po c -lu)t 


(17) 

(18) 


then 

tp° x + = o (is) 

Thus, the governing equation for the velocity po- 
tential becomes Laplace’s equation, which indicates 
a diffusive propagation rather than a wove like 
propagation. The effect of this propagation on p 
and the acoustic velocities will now be considered. 


The continuity equation, Eq. (13), simplifies 
by the use of Eq. (17) to 


ft + U ^“° 


( 20 ) 


Again, nsaumlng that p « p 0 0 -ia*, £q, (20) yields 


pO+i I¥° 0 


(21) 


For a plane pressure wave at x «* 0, a solution for 
p° from Eq , (2i) Is b imply 

p o „ Cie iarc/t; (22) 

Similayily aaauming u » u°e’^ u1 *' and v ® v £> o'‘^ ,U3 *' 
and substituting Eq. (22) Into Eqs. (14) and (15) 
yielda after inegration 
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These solutions of the continuity and momentum 
eq, ations have lead to the fact that 


p, U ,V,aeCi<Vu)x (25) 

that is, only forward propagating perturbation 
pressure p and velocities will exist and they 
will propagate at the speed of the medium flow U. 
Therefore, Eq. (16) does not have an acoustic solu- 
tion, consequently, the potential flow solution 
cannot be used to estimate the effects of shear 
even in an approximate manner. 

For the case of boundary layer flow, however 
a procedure based on the Golds tcin-Ricc anatysisl4 
might be used for predicting the attenuation with 
shear flow in conjunction with the potential finite 
element method for turbofnn inlets. For a linear 
sheer layer, the impedance at the wall can be 
transformed into an effective impedance, Z ej -£, flt 
the edge of the boundary layer. Tills effective 
impedance could be used as the boundary condition 
for the potential flow calculation. The procedure 
is approximate because the calculation assumes a 
uniform boundary layer in the axial direction. 

For n growing boundary layer, the transform could 
be assumed to apply locally. In any event, this 
procedure represents a simple affective way of ac- 
counting for shear with a potential flow program. 


Attenuation Sensitivity to Variable Area 

The theory used for the design of acoustic 
liners for turbafan inlets without splitter rings 
has lined a cylindrical duct for a model. As men- 
tioned earlier, the present paper will estimate the 
sensitivity of duct attenuation to inlet curvature, 
centerbody and flow gradients (irrotationul) for a 
typical turbofan inlet. To accomplish this task, 
a number of acoustic calcu Lotions are performed on 
tile Boft-walled turbofnn inlet shown in Fig. 1. 
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Thio in lot geometry is oirallor to the QCSEE (quiet, 
clean, short-haul experimental engine) inlot which 
was designed to suppress inlet- emit ted engine ma- 
chinery noise using a high throat Mach number mid 
thus has somewhat mora curvature than a ’.unvention- 
ui inlet, in this esao, the velocity potential 
finite element program of Kel. 1 was taudifiod to 
include quadratic elements. 

The governing equations in the program are de- 
veloped from the axlsyinaetric cylindrical coord- 
inate form of Eqs. (1) to (6), The exact form of 
<P 1b 

9(r,0,i,t) - (26) 

The detailed govarnlng equations will not be pre- 
sented here since they are fully documented in 
Ref. 1. Boundary layers are not considered in the 
calculation because, as was just shown, the program 
is limited to only lrrotstlonal flow fields. 

At the present time, tiio frequency range and 
radial mode content Is limited by the maximum num- 
ber of elements which tho computer con handle. The 
analysis considers plane wave input with dlmcnolon- 
leau frequencies of q up Co 2. In tills case, 
seven elements wore used to resolve the radial 
modes and fourteen elements wore UBCd to resolve 
axial variations. The acoustic design of an actual 
inlot will, of course, consider much higher fre- 
quencies. Unfortunately, at the present time finite 
element schemes can not handle high frequencies , 
Therefore, this analysis only investigates low fre- 
quency noise in variable area inlets. This might he 
corrected by adapting tho wave envelope technique of 
Ref, 6 to this problem. The wave envelope technique 
can handle high frequencies with much fewer elements. 

Attenuation calculations for a variable area, 
inlet, a straight open cylinder and a bellmouth 
(Fig. 2) all with on L/l) of l nre presented In 
Fig. 3. The ratio of treatment length to duct 
diimeter in these cases is 0.5. The soft wall ox- 
t nds over half the length of the duct beginning at 
the nose of the centorbody, With this guometry, 
the infinite uniform cylindrical duct model could 
be expected to give u reasonable estimate of the 
attenuation. The results displayed In Fig, 3 ore 
for no mean flow. The effects of flow will be con- 
sidered later. The impedance values chosen ere the 
some for each inlet at a particular frequency. 

They are listed in Fig. 3, The impedance values 
chosen in the calculation approximate the optimum 
Impedance (max. attenuation) for a plane pressure 
wave input into a long circular duct. A plane 
acoustic velocity is used nu the source with a 
p o c 0 exit impedance. 

The effect of inlet curvature and the center- 
body on the duct attenuation cun be ascertained by 
comparison of the various curves in Fig. 3 repre- 
senting each inlet. The cylinder and bellmouth in- 
lets give approximately the same attenuation as o 
function of frequency. Since the bellmouth is of 
uniform circular geometry, the uffoct of the short 
centerbody on the results appears to be insignifi- 
cant in this no flow situation. The large area 
variation of the other inlet, however, does signif- 
icantly affect the results. Tills suggests that re- 
flections in this inlet may be important. 

Figure 4 compares preliminary attenuation cal- 
culations for the vorlable area inlet and a cylin- 


drical inlot in a flow situation. Hie fan plane 
Mach number for tho variable area inlot was 0.52 
with an average Mach number of 0.579 ot the aofb 
wall dcction of the inlet, as calculated by ane- 
dlmcnsional iscntropic gas dynamic relations. The 
Msch number in tho cylinder was also choson to bo 
0.579. The impodanco values used in these calcula- 
tions were these listed in Fig. 3 divided by 
(1 - M) 2 with M cquul to 0,579. Aguiti, us seen 
in Fig. 4, the QCSEK inlet differs considerably 
from the cylinder results. 

In Fig. 4, the flattening at the cylinder at- 
tenuation curve probably results because tho Imped- 
ance values chosen are not at the optimum for this 
Mach number, Howevur, the important point la that 
for the some value of impedance, area variations 
can significantly change the calculated attenuation, 


Far Field Radiation 

As mentioned in the introduction, this section 
will be concerned with current ideas for extending 
finite elements from the internal portion of the 
duct into the far field. Dndur a three year exten- 
sion of NASA Lewis grant NAS-3036, which will run 
from 11/7B to 10/Si, profeoaors 2inn and Sigman 
will be developing metllodB for extending their fi- 
nite element onalyula into the far field. First, 
a small effort will be devoted to more accurately 
and efficiently specifying the steady flow field 
around an inlet in a form compotnbie with the acous- 
tic calculations, However, the majority of the 
work effort in this grant will bo devoted towards 
the calculation of «countic properties of the inlet 
and the associated external region. 

The direct extension of finite elements into 
the fur field is illustrated In Fig. 5. The num- 
ber of elements greatly Increases the further the 
analysis is carried from the inlet. Therefore, 

Ibis approach requires very efficient woys of stor- 
ing and solving the finite element global matrix. 
Conaiderable effort should be devoted to the direct 
method, since improved handling of the elements 
will benefit far field approachea as wall os im- 
prove canobilltios for present Induct analysis. 

As mentioned earlier, to handle higher frequencies 
inaide the duct, greater efficiency is still needed 
to increase the number of internal elements. 

A new numerical partition method partitioned 
matrix approach is now presented which allows far 
field acoustic problems to be subdivided into 
smaller Independent and thus more raanngable prob- 
lems. This procedure can be incorporated in previ- 
ously developed programs with minimal effort. Sam- 
ple calculations for a two-dimensional rectangular 
coordinate soft wall duct with uniform flow are 
presented later to illustrate the method. Before 
discussing the partition approach, the boundary 
conditions used in a typical numerical acoustic 
analysis will be briefly reviewed. 

Finite element or finite difference analyses 
usually require an acoustic pressure or velocity 
distribution us thu entrance boundary condition and 
some assumed impedance as the exit boundary condi- 
tion. For plana wave propagation without mean flow 
in a finite Btraight duct, the acoustic impedance 
at the outlet of the duct is nearly p 0 c 0 for 
ditnensicmiesB frequencies q greater than 1.5, as 
shown in Ref. 16. Figure 6(u) depicts the geometry 

.iginal page is 
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and appropriate boundary conditions for a rectangu- 
lar duct with these conditions. For spinning or 
higher order modes, the exit Impedance values of 
Refs. 1 or 3 could be used. 

The basic idea for the partitioning approach 
comes from Che marching technique developed in 
Ref. 17. As shown in Ref. 17, in regions where re- 
flections are small, the P 0 c a exit impedance was 
nearly constant along the entire length of : le duct. 
Thus, the internal portion of the duct and tile far 
field can be separated Into regions ond solved in- 
dependently, The exit pressure calculated in the 
first region is used as the initial condition for 
the pressure In the second region, and so forth. 

For high frequency sound where reflections are 
small, the partition could bu taken at the exit, aa 
shown ill Fig. 6(a), For low frequency sound or for 
an arbitrary input with multiple nodes, the parti- 
tion should be moved from the exit to the far field 
whore the f> 0 c 0 Impedance is valid, as shown in 
Fig. 6(b). 

The number of partitions used in the far field 
could be as many as desired, as illustrated in 
Fig. 7, The total computer storage is reduced by 
approximately the square of the number of parti- 
tions, If seven partitions were used, this would 
be a storage reduction of about 50. Furthermore, 
since the solution times ure roughly proportional 
to the total number of nodes cubed (Ref. IS, 
pg. 261), the total running time should also bu re- 
duced by a factor of 50. 


To illustrate the capabilities of the parti- 
tioning technique, the noise attenuation at the 
optimum Impedance (point of maximum attenuation in 
the Impedance plane) is calculated using finite 
difference theory for a two dimensional rectangular 
duct with an L/li of 3.63, a uniform Mach number of 

0.3 and a plane wave input. 

Dimensionless frequencies of jj equal to 1.5 
and 2 ure considered. As shown in Ref. 16, a di- 
mensionless frequency of 1.5 Is Just sufficiently 
high so that a p 0 c o exit impedance is valid. In 
these two cases, ns discussed in Ref. 17, the p 0 c 0 
exit impedance is nearly constant along the entire 
duct. Therefore, the duct can be partitioned. 
Figure 8 allows the duct with seven assumed parti- 
tions represented by the dush vertical linos. 

Figure 9 displays the calculated maximum sound 
power attenuations for q equal to 1.5 and 2. The 
number of axial grid points K have been varied to 
check for convergence. The results converge to the 
analytically predicted attenuation 1 ^ which apply to 
ducts of infinite length. Without partitions, ap- 
proximately 100 nxial grid points are required 
while with partitioning only 20 axial grid points 
are needed in each partitioned element. As men- 
tioned in the body of the report, this reduces the 
storage requirement by a factor of 50 and the total 
running time also by a factor of 50, Of course, 
the method must be tested over a greater range of 
acoustic variables to check its validity. 


Concluding Remarks 

The finite element velocity potential program 
developed for NASA Lewis by Georgia Institute of 
Technology Is suitable for handling low frequency 
acoustic wave propagation in ducts with area varia- 
tions, centerbodles, and axial variations in wall 
Impedance whenever irrotatlonal flow gradients 
exist in a duct, Consequently, the program Is 
ideally suited for handling low frequency acoustic 
wave propagation in a turbofan inlet in which the 
moan flow field la predominantly irrotatlonal and 
two dimensional, 

Unfortunately, us shown herein, the velocity 
potential program can not bo used directly to esti- 
mate the effects of wall shear layers on acoustic 
propagation, However, a boundary layer correction 
based on the Goldstcln-Rlcc analysis is suggested 
for predicting attenuation in soft wall ducts. 

In some sample calculations for a variable 
area inlet, the combined effect of inlet curvature, 
and flow gradients was shown to have significant 
effects on the attenuation of a given acoustic 
liner for very low frequencies. The short center 
body did not seem to have such effect. 

Finally, approaches are discussed that are 
actively being considered for extending the finite 
element solution to include the far field os well 
as the internal portion of the duct, A new matrix 
partitioning approach is presented that can be in- 
corporated in previously developed programs which 
may allow the finite element calculation to be 
marched into the for field. The partitioning ap- 
proach provides o large reduction in computer stor- 
age and running times. 
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Figure 2. - Bellmouth inlet. 
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Figure 3. - Sound attenuation for inlets with a treatment 
UD 0 of 0. 5 for a plane acoustic velocity source with 
zero mean flow. 
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Figure 4. * Sound attenuation for inlets with a treatment 
L/D 0 of 0. 5 for a plane acoustic pressure source with a 
fan face Mach number of 0. 52. 



Figure 5. * Extending finite element into far field. 


P 0 > r REGION l_ D c REGION NO. ? 

X* >| p °° , 


(a) HIGH FREQUENCY APPROXIMATION. 
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(b) LOW FREQUENCY APPROXIMATION. 

Figure 6. * Application of partitioning to engine inlets. 
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Figure 7. * Multi -partitions used in tar-field analysis. 
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Figure 8. - Two-dimensional duct partitions. 
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Figure 9. - Effect of number of axial grid 
point K on attenuation for seven 
partitions in a soft wall duct UH • 3.43, 
M • 0. 3 and wall impedance at optimum 
value for plane wave input. 


